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TECHNICAL NOTE 2658 


LAMINAR BOUNDARY LAYER OVER FLAT PLATE IN A FLOW HAVING 

CIRCULAR STREAMLINES 
By Artur Mager and Arthur G. Hansen 


SUMMARY 

The laminar-houndary-layer development over a semi-infinite flat 
plate placed in a flow with .concentric circular streamlines was investi- 
gated "with the limitation of small total turning of the main-stream 
flow. The shape of the velocity profiles in the direction of the main- 
stream flow and perpendicular to it was analytically determined for an 
incompressible flow and a compressible flow with Prandtl number egual 
to 1. 


The boundary- layer thickness was shown to be proportional to the 
square root of the distance from the leading edge of the plate when 
measured along the streamline of the main- stream flow. The deflection 
of the boundary-layer flow at the plate surface from the direction of 
a circular streamline in the main flow was shown to vary directly with 
the turning. With increase in the Mach number of the main-stream flow, 
both the boundary-layer thickness and the deflection increased. 


INTRODUCTION 

The relative lack of theories explaining the behavior of boundary 
layer when a lateral curvature of the main-stream flow exists has been 
especially apparent in the application of aerodynamics to the design of 
turbomachinery. The development of the boundary layer in such Gases is 
strongly influenced by the corresponding normal pressure gradient toward 
the center of curvature, giving rise to a component of "secondary flow" 
in the boundary layer. For the laminar case, most of the published work 
has been restricted to yawed cylinders, wings, and cones (references 1 
to 4). One notable exception of direct application to the design of 
compressors and turbines is reference 5, wherein the boundary layer on a 
rotating blade is analyzed. For the turbulent case, a general but 
approximate solution of the momentum- integral equations based on an 
assumed velocity distribution and friction law is obtained in refer- 
ence 6. 
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Aside from the conventional boundary-layer approach to this prob- 
lem, a number of investigators (references 7 to 9) have obtained solu- 
tions for secondary flow arising from flows of varying total pressure 
or varying enthalpy by neglecting the influence of viscosity but admit- 
ting the existence of vorticity. Although such procedure obviously does 
not permit satisfaction of all the boundary conditions (because the 
order of the general differential equations for the flow is reduced), 
the results so obtained give a fair check with the experimental data, 
except in the regions close to the wall. Thus the indications are that, 
while it is possible to obtain a fair picture of three-dimensional flows 
in the preceding cases by neglecting viscosity, such procedures are 
inadmis sable where thin boundary layers exist. 

The object 'of this investigation, which was conducted at the TTACA 
Lewis laboratory, is to solve a case somewhat analogous to that of 
references 7 to 9 in regard to main -flow orientation while retaining the 
usual Prandtl boundary-layer assumptions, and thus to demonstrate the 
inf luence of viscosity and compressibility on the secondary flow in 
relatively thin boundary layers. 

To this end, the main flow outside the boundary layer is assumed to 
follow concentric circular streamlines in planes parallel to a semi- 
infinite flat plate and to be uniform in a direction perpendicular to 
the plate. The edge of the plate is placed so as to coincide with a 
radial line through the axis of the main flow (fig. l). It is assumed, 
of course, that in the established main-stream flow the effect of vis- 
cosity is neglibible. In the simplification of the equations for the 
flow, it was found convenient to follow the procedure first used by 
Moore (reference 4) . 

Although the solution is primarily of theoretical interest, it may 
have some direct bearing on the understanding of secondary flows over a 
portion of the hub between the compressor guide vanes where the boundary 
layer is probably still thin and laminar, and where the main flow under- 
goes relatively little turning, often along a circula r path. Here again 
the analysis would probably be applicable only in regions where boundary- 
layer flows from the blade surfaces and other "wall effects" were not 
influential. 
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BASIC EQUATIONS 

The following equations describe the motion of a compressible vis- 
cous fluid in vector form (reference 10): 

Conservation of mass 


|+pN=o (i) 

Conservation of momentum 

P It * P* “ VP + MV 2 ! + | MV (v-l) + 2 J"(VM) -v] 1 + (vt-Oxfcxq) - 

| (v!)(vm) (2) 


Conservation of energy 


P §f + p(v*g) = A + V- [(kV)T] 

where the dissipation function is given "by 

A = M jsv- [(i*V)qJ + (VXq) 2 - 2q-v(?-q) - | C^-a) 2 ^ 


Equation of state 


(3) 


p = PRT (4) 

(All symbols are defined in the appendix.) 

When dealing with problems involving a lateral curvature of the 
potential flow it is convenient to use an orthogonal coor dinat e system 
illustrated in figure 2. This system is characterized by a fixed refer- 
ence axis x of arbitrary curvature c = c(x) and is related to a 
Cartesian system as follows: 
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Xi 



cos p dx + z sin (3 


A 


x 2 = y 


( 5 ) 


Xg = constant + z cos p - / sin p dx 

JO 


J 


where 



P - P(x) 

The elements of length at (x,y,z) in the direction of the increas- 
ing coordinates are 


hpdx h 2 dy hgdz 

and their values may he obtained by use of equations (5) as 


hj = (l + cz) 


h 2 = ! 
h 3 = 1 

The transformation of equations (l) , (2), and (3) from the Cartesian 
system to the orthogonal system is readily obtainable by use of the 
expressions for the elements of length (see., for instance, reference 11 
or 12) . When the y-axis is assumed perpendicular to the surface over 
which the flow takes place and the usual Prandtl assumptions about the ' 
order of magnitude of various terms are made, the equations (l), (2), 
and (3) take the following form for steady flow in a boundary layer where 
equations (6), ( 7 ), and (8) correspond to equations (l), (2), and (3), 
respectively: 


2357 



2357 


NACA OUT 2658 


5 


[iz ( pu ) + [c 1 + cz ) pv ] + h t 1 + = 0 


(6) 


/ u 5u Su du c 

p \^1 + cz 5x + v 5y + w §z" + 1 + cz uw 


r ) = -^ 


1 9p d ( Bu \ 
cz 5x + 5y §y/ 


(7a) 


( u dw , 3w , 

\przz s + v + 


„ dw c 2 l 3p , d / cftr\ s 

w §7 " TT^£ u J = " + ^ V p ’Sj) ( n) 


Concomitant with (7a) and (7h) is the result that ^ is of order of 
magnitude 5 ~ 0. 

( u cSD dE cffi\ p Jchi ^ r/_ . i 

p \ — ™ 5i + v ^ + w 5iy = - r+^\£E + K 1 + cz H + 

^ [(i + cz)w]| + +(|p) + |f ( k §|) (8) 

The equation of state (4) remains unchanged. 


Because of the nature of the problem under consideration, several 
simplifications in the preceding equations can he made. First, the fixed 
reference axis x will he chosen to coincide with a streamline of the 
main-stream flow where these streamlines are assumed to he arcs of con- 
centric circles in planes parallel to a semi- infi nite flat plate (see 
fig. l). This flat plate is oriented so that its leading edge is orthog- 
onal to the reference axis xj hence the leading edge can he defined by 
x = 0, y = 0. Secondly, if the velocity of the fluid in the main stream 
is designated uq, it follows that uq = constant along any given stream- 
line; that is, uq = uq(z) . 


It also follows from the nature of the prescribed main-stream flow 
that the pressure gradient in the x-direction and outside the boundary 
layer is identically zero. In view of the order of magnitude of Bp/By, 
the same pressure prevails In the boundary layer as in the main flow, 
and consequently (where the subscript indicates partial differentiation) 


p x = 0 


( 9 ) 
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In addition, for flow with, circular streamlines it may "be shown 

that 


~ ( 1 + cz ) 

The bo undar y conditions on equations (6) and ( 7 ) are 


( 10 ) 


u(x,ccr,z) 
W(x,00,z) 
u(x, 0 , z) 

where 

and at the leading edge 


= u x (z) 

= 0 

= w(x, 0,z) 
x > 0 


A 


v(x,0,z) = 0 


> 


u(0,y,z) = u x 
w(0,y,z) = v(0,y,z) = 0 

where 

y > 0 > 


(H) 


For equation (8), usually two boundary conditions on temperature are 
used by def ining its value at y = 0 and y a>. For the case under 
consideration, however, it will be necessary to define temperature only 
at y = 0. 


SOLUTION OF EQUATIONS 

Application of Vector Potential Concept 

In order to satisfy the equation of continuity (6) and to eliminate 
one dependent variable from equations ( 7 ), the concept of a vector poten- 
tial is employed (see reference 4 ) . To this end, two functions t)r and 
9 are defined, where 
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Pu=* y 

(l + Cz)pW SC Py 

(1 + ez)pv s -(ijr x + cp z ) 


> 

J 


( 12 ) 


It may "be verified "by direct substitution of equations (12) into equa- 
tion (6) that equation (6) is identically satisfied. 


Substitution of equations (9), (10), and (12) into equations (7a) 
and (7b) yields 


*y ♦ 


^ + (v)y + r hz ty Q 


(l + cz) 


♦y 

LI l — 

V P/7 


(13a) 


% 


" qy 

p(l + cz) 


X 


<p 


y 


(l + cz) 


y 


- ^ 


p(l + cz) 


Jy 


f %) 2 


(l + cz)^ P 


qy 2 

2 = cp l U l + 


P 


Sr\“ 

p /y 


(13b) 


Incompressible Case 

The solution for the incompressible flow (p = constant) will be 
presented first. For this case it is necessary to consider only equa- 
tions (13) in the boundary layer; the energy equation need not be used. 

In accordance with the usual procedure in the solution of incom- 
pressible boundary-layer flows, the following substitutions for the 
variables are employed: 


/ 
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/ U 1 ^ 
^ ^ / Vx(l + zc) 

ijr = P [/xUu 1 (l + zc) F(t]) 

Cp = pcx (l + zc) i/xVuq^I + zc) G(t]) y 


(14) 


The following system of equations is obtained by substituting 
expressions (14) along with their appropriate partial derivatives into 
equations (13) (note that p/p = V ) : 


(cx) 2 G’F 


(l + zc) 

U-^C 


dU]_ 

dz 



/ 


(ox) 2 F"G 



(1 + zc) ^1 
2cuj dz 


FF” 

~ 


- F’ ' 1 = 0 


(l5a) 


F'G* - ^ + (l - F' 2 ) - G"’ + (cx) 2 G' 2 ~ 


(cx) 2 G"G 



(1 + zc) du ll 0 
2cu^ dz J 


It is interesting to note that setting 

U 1 = ^oC 1 + zc ) n 


(15b) 


(16) 


in (15a) and (15b) reduces these equations, respectively, to the forms 


(cx) 2 G*F* (n + 1) + (cx) 2 F”G 


! +n 


FF” 

~Z~ 


_ F* ' ' = 0 


(17a) 


n(cx) 2 (G*) 2 - (cx) 2 G"G^§ + nj + F*G' - + (l - F' 2 ) - G'." = 0 

(17b) 
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Letting n = -1 in equation (16) gives the expression for the velocity 
in the classical irrotational vortex flow, whereas lett ing n = 1 
yields the expression for a special case of "wheel-type" flow. 

How F and G and their derivatives are assumed of order of mag- 
nitude of unity, while (cx) 2 « 1. It is seen that cx is the total 
turning of the flow outside the houndary layer, and because of the fore- 
going restriction it must remain small. It is further assumed that 

du]\ 

1 is of the order of magnitude of uni ty (which should he the 

case for most flows of interest and is evidenced in part hy equation (16) ) . 
Heglecting all terms of the order of (cx) 2 in equations (15) and (17) 
reduces these equations to the following two total differential equa- 
tions: 


( 1 + zc 
cui 


FF" + 21" * * = 0 (18a) 

2F‘G* - FG" - 2G»«* + 2[l - (f 0 2 J 
The corresponding boundary conditions are 

F = F* = G = G* = 0 at T] s 0 

and 

F* = 1 and G' = 0 as r) — » oo 

Equations (18) and (19) show that for small turning F is the 
usual Blasius function and that equation (l8b) may he solved for G' 

(or G) since F is known. Such solution was obtained hy relaxation 
procedure i n asmuch as F and its derivatives are given in tabular form 
only (cf. reference 13, p. 121). The results from this computation are 
presented in figure 3 and table I. 


= 0 (l8b) 


> (19) 

J 


Boundary-layer thickness . - If the boundary-layer thickness 5 is 
defined as the value of y at which u is a certain percentage of u^ 
(usually over 99 percent) , a corresponding value for T] may be found 
from the solution for F where this value of tj is designated rj*. 
Consequently the expression for boundary-layer thickness may be obtained 
from equations (14) as 


5 = 



( 20 ) 
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or nondimensionally in terms of the reference length x as 

_1 

| = n* + cz Ee 2 (21) 

On a given streamline, it can "be seen from equation (20) that the 
boundary-layer thickness varies directly -with the square root of x. 
For the special case of free vortex flow, equation (20) becomes 

5 = T]* (1 + cz) (22) 


and. it follows that for any given x and uq^ thickness increases in 
proportion to (l + ez) or that the boundary layer becomes thinner with 
an increase in curvature of the flow. The lines of constant thickness 
are obtained, by considering 5 to be a constant in equation (22), or 

x = ( 23 ) 
(l + cz)^ 

Therefore, the curves of constant 5 are spirals of JLituus; they are 
shown in figure 4. For more general flows, the lines of constant 
boundary-layer thickness are given by 

u-^ (constant) 
x - (1 + zc) 

where, of course, uq is a function of (l + zc) . 

Flow direction at plate surface . - From equations (12) and. (14) the 
egressions for u and. w are 


U = Uq F* 


w = cxujG' 


Therefore, if 


e 


w 

iim ~ 
T]-* 0 
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the flow direction at the plate is given by 


a = arc tan e 

However, since G* (o) = F'(o) = 0, 


r G"(o) 

cx fW 


Equation (25) indicates that 


inasmuch as 


T — e T 
0,z - 0,x 


T _ pBu 

°’ x ~ Srj(y = 0 ) 


Bw 


T °’ z & J(y = 0) 

The ratio G"(0 )/f"( 0) was evaluated numerically to give the result 


( 24 ) 


(25) 


€ — -3.26 cx 

It may "be seen, then, that at the plate surface the flow is deflec- 
ted from the direction of the streamline toward the center of curvature 
of the main-stream flow and that this deflection is proportional to the 
turning. The lines of constant e are seen from equation (25) to "be 
those of constant turning, that is, radial lines (see fig. 4). 


Effect of Compressibility 

Complexity prevents the extension of the previous results to include 
the effects of compressibility in the same manner as for the incompres- 
sible case. For compressible flow, a solution can "be obtained only in 
the immediate neighborhood of a reference stre amlin e on which the Mach 
number has been specified. However, it is possible to calculate the 
Mach numbers on various chosen streamlines in terms of the Mach number 
specified on an initial streamline. Solutions can then be found in the 
neighborhood of the chosen streamlines using the calculated Mach numbers. 
To find the boundary- layer flow in any region of interest is therefore 
possible . 
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Simplification of energy equation . - Letting 6 - C p T and assuming 
that the specific heats and Prandtl number are constant, the energy 
equation can be written in terms of cp and t}t as 


V 


+ cp 0 

x y z 


- Ct 


x + *z> 0 y 


c 

1 + cz 




+ (l + cz)p 



ft 




p(l + cz)J 



(26) 


For Prandtl number equal to 1 and zero heat transfer at the plate sur- 
face, equation (26) is replaced by 


CpT + 



(27) 


where T w is the constant temperature of the plate. 

Keduction of equations to nondimensional form . - The following 
approximate formula (reference 14) is used to represent the relation 
between the viscosity and the temperature: 


-&--A 

“ 1,0 


T 


T 


1,0 


(28) 


where A is a constant suitably adjusted to give the best agreement 
with the actual relation over the temperature range considered. 


It now becomes convenient to make all physical quantities in the 
various equations nondimensional. The following system suggested in 
reference 4 will be used; all quantities on the left are to be con- 
sidered relative to the quantities on the right, which are measured 
outside the boundary layer and along the reference axis x: 


u,v,w 

relative 

x,y,z, | 

relative 

p 

relative 

T 

relative 


to u^o 

to pi,oVpi,o u i,o 

to p^o 

to u2 i,o/ 2C p 
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p relative to Pl,o u2 l,0 

■'I'lCp relative to ii^qA 

All subsequent equations will appear in nondimensional form. As a 
consequence of the use of the nondimensional form, equation (13a) 
becomes J 


t 




y \ p 


h 

p 


Ctx + v 


Cf): 


y ‘ 1 + zc ^y 



= (l + cz) 


T 


T 


1,0 


th) 

\ p / y. 


equation (l3b) becomes 





yy 

(l + zc) 



- (t- 


x + < P z ) 


0 


<Py 


p(l + cz) 


£ Air c 


T (O - 7 -p = -Cp-.Un*- + 

p y (l + cz) 2 P 1 1 


2 O. f E ^ " 
_ T 1,0 VP /y. 


equation (27) becomes 


T + < 


©’•biftj}- 


5V 


and equation (l) will be written 


Ai- = p _ 2 L_ 

* 1,0 * 1,0 


( 29 ) 


Transformation of Howarth . - The following transformation of vari- 
abies indicated by Hcwafth (reference 15) is now utilized to mate the 
nondimensional form of equations (13a) and (13b) s imilar to the incom- 
pressible form of these equations: 
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Y = 



pdy 




X = x 


z 

* 


z 


Pi 

* 1,0 


¥ 


cp S / P l 


1 , 0 . 


$ 


> 


J 


(30) 


As a consequence of transformations (30), equation (29), and the 
fact that Bp/dy is of the order of magnitude 5, equations (l3a), 
(13b), and (27) become, respectively: 


¥ Y ¥ YX + $ Y ¥ YZ - O^x + $ z) ¥ YY + 1 + cz ¥ y $ y 
= ¥ yy (rfcS ^ ^ i^*) + (1 + cz)¥ ^ ' (31a) 

¥ Y%X + $ Y $ YZ - (^x + $ Z^YY " l + cz ^Y^ + C ( 1 + cZ ^ [fp - U 1 " ] 


c ^l 2 T l,0 S> * f 
rr^Z ~ ^ i^YY 1 CZ;W YYY 


(31b) 


T + 



(32) 


The same procedure as for the incompressible flow case is now 
followed. The similarity variable and corresponding vector potential 
components are defined as 
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T) = Y 




U-, 


X(l + Zc) 


¥ = /xu-l( 1 + Zc) F(t]) 


® - cX(l + Zc) ^Xuj_(l + Zc) G(t|) 
By use of these definitions , equations (3l) hecome 




(CX) 2 G.?. [l + 5i] + (oX)2s r„ j-| 

(1 + Zc) du l 1 

il dZ J - 


rru»Jt I IT lx 2 3 


M- 


L 1 . 


dU T *1 FF" _ 


2cUn 


- F* » * =0 


(33) 


(34a) 


P'G' + + W ! . f,2 1. o„, + 


(cX) ^G 


2 P 1 2 


(i + Zc) du i r-1 

L cutl dZ " 2 M 1 


- (cx) 2 G"G 


1 + ( X + Zc ) + I M 2~ 
.2 + 2cun dZ + 2 M 1 


= 0 


(34b) 


Boundary conditions are the same as in equations (19) with the 
appropriate change in variables) that is 


F=F' = G = G' = 0 at tj = 0 


and 


F 1 = 1 and G* = 0 as T) — 5 > co 

Once again when terms of the order of magnitude of (cx) 2 are 
eliminated under the same type of restrictions as in the incompressible 
case,, equation (34a) becomes the Blasius equation for F. However, it 
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may lie seen that equation (34b) cannot; "be solved directly as a total 
differential equation in rf since Mq is a function of Z. In fact, 
if uq is specific ally prescribed as a function of (l + cz) (that is, 
Uq = f(l + cz) where f = 1 when z = 0), then Mq may he expressed as 



anfl in the immediate neighborhood of the reference axis X, Mq « Mq^ q 
and may he considered a constant. Incorporating the preceding restric- 
tions into equation (34h) gives 

F'G' _ g«m _ q/ 2 FG" + (l + jl - (F') 2 J = 0 (36) 

This equation has the solution 

G' = ^1 + 2 ^1, 0 J C 1 !) (37 ) 

where G 1 is known from the incompressible case. Curves of G 1 (tj) for 
various values of Mq^g ajre plotted in figure 5. 

As mentioned previously, the flow in the neighborhood of stream- 
lines other than an initial reference streamline may he calculated in 
the same manner hy taking the reference axis along the new stre aml i n e 
and hy computing the appropriate Mach number on this streamline from 
equation (35) . 

Physical coordinates . - Relations (29)_, (30), and (32) may he used 
to obtain values of y as a function of tj for a fixed value of x 



It is to he recalled, of course, that these coordinates are the non- 
dimensional form of the true physical coordinates. 

Interpretation of results . - From equations (18) , (36), and (38), 
the effect of compressibility may he seen to he twofold, namely, changes 
occur in the coor din ate y as a function of T) and the differential 
equation for G. These changes have a corresponding influence on the 
velocity distribution, the houndary-layer thickness, and the flow deflec- 
tion at the plate surface. 
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Since 


n 


s: 


(f*) 2 an 


is a positive function of t), for any 

_z “ _i 

given value of rj the value of y increases with Mach n umb er, as shown 
hy equation (38). Consequently, the houndary-layer thickness 5 
increases with Mach number. It also follows for a chosen Mach number 
that 


du 



r-i 

2 





thereby indicating that du/dy increases with Mach number for 
tj = constant. The same effect is known to exist in two-dimensional 
boundary- layer flow. 

As far as the flow deflection at the plate surface is concerned, 
the changes in the shape of G' (see fig. 5) and equation (25) indicate 
that for a given turning the deflection toward the center of the flow 
field increases with Mach number. Consequently, larger secondary flows 
in the boundary layer are obtained with the increase in Mach number. 


CONCLUSIONS 

The following conclusions were drawn from a theoretical investiga- 
tion of boundary-layer flow over a flat plate when the main- stream flow 
outside the boundary layer follows concentric circular streamlines and 
the tota.1 turning of that flow is small: 

1. The profile of the velocity component in the direction of the 
flow .outside the boundary layer is identical to that existing in two- 
dimensional flow over a flat plate. 

2. The magnitude of the velocity component in the direction perpen- 
dicular to the flow outside the boundary layer and in a plane parallel 
to the plate varies directly with the turning of the main- stream flow. 

3. The boundary-layer thickness varies with the square root of the 
distance along the streamline of the main-stream flow as measured from 
the leading edge of the plate. For the special case of free vortex 
flow, the lines of constant boundary- layer thickness are spirals of 
lituus. 


4. The deflection of the flow at the plate surface from the direc- 
tion of the main-stream flow is directly proportional to the turning of 
the main-stream flow. The lines of constant deflection are therefore 
radial. 
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5. Both the Boundary-layer thickness and the deflection at the 
plate surface increase with Mach number. Consequently, larger secondary 
flows in the boundary layer are obtained for higher Mach numbers. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November 27, 1951 
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APPENDIX - SYMBOLS 


The following symbols are used in this report: 

A constant in viscosity-temperature relation 
Cp specific heat at constant pressure process 
c curvature of x-axis, df3/dx 


JD_ 

Dt 

V 

E 


v,v\ 

G,Gj 


Eulerian derivative 

vector differential operator 
energy 

functions related to components of vector potenti al 


f inertial field force 

h-i factor used in transformation of coordinates (i = 1, 2, 3) 

k thermal conductivity 

M Mach number 

p pressure 

q. velocity vector 

R gas constant 

Re Reynolds number hased on x, (u^x/v) 

T temperature 

u,v,-w velocity components in curvilinear coor dinat e system 

H,T,Z coordinates defined "by equation (30) 

x,y,z curvilinear coordinates 

a "boundary-layer deflection angle measured from direction of 

flow at plate surface to direction of flow outside 
"boundary layer 



20 


UACA BJ 2658 


P 

r 

8 

e 


ii* 

e 

A 


p 

v 


p 

T 0,x, T 0,z 



angle between X^-axis and tangent to x-axis 
ratio of specific heats 
houndary-layer thickness 

measure of houndary-layer deflection near plate surface, 
tan a 

similarity variables 
value of T] at y = 8 
enthalpy, CpT 
dissipation function 
absolute viscosity 
kinematic viscosity, p/ p 
density 

shear stress at 'wall in x- and z-directions, respectively 
components of vector potential 




Cartesian coordinates (i = 1, 2, 3) 


Subscripts: 


1 outside the boundary layer 

0 on reference axis x 

w at plate surface 

All other subscripts indicate partial differentiation. 
Superscript primes indicate total differentiation. 
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Figure 1. - Flow 
orientation of 
layer. 




Figure 2. - Transformation from Cartesian coordinates 
to curvilinear coordinates x,y,z. 
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Figure 3. - Functions F’ and -G-' plotted against 
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Figure 5. - Effect of Mach number on function G'. 
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